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Core-collapse supernovae emit of order 10^* neutrinos and antineutrinos of all flavors over several 
seconds, with average energies of 10-25 MeV. In the Sudbury Neutrino Observatory (SNO), which 
begins operations this year, neutrinos and antineutrinos of all flavors can be detected by reactions 
which break up the deuteron. For a future Galactic supernova at a distance of 10 kpc, several hundred 
events will be observed in SNO. The and Vr neutrinos and antineutrinos are of particular interest, 
as a test of the supernova mechanism. In addition, it is possible to measure or limit their masses by 
their delay (determined from neutral-current events) relative to the Ve neutrinos (determined from 
charged-current events). Numerical results are presented for such a future supernova as seen in SNO. 
Under reasonable assumptions, and in the presence of the expected counting statistics, a Vfj, or Vt 
mass down to about 30 eV can be simply and robustly determined. If zero delay is measured, then 
the mass limit is independent of the distance D. At present, this seems to be the best possibility 
for direct determination of a i/^ or Vr mass within the cosmologically interesting range. We also 
show how to separately study the supernova and neutrino physics, and how changes in the assumed 
supernova parameters would affect the mass sensitivity. 

14.60.Pq, 97.60.Bw, 25.30.Pt, 95.55.Vj 
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I. INTRODUCTION 

As emphasized by Weinberg |]l| and many others, 
whether or not neutrinos have masses or other properties 
beyond the standard model are questions which address 
some of the deepest issues in particle physics. Yet almost 
seventy years after they were proposed by PauU, most of 
the properties of neutrinos are defined only by limits. 
In particular their masses, if any, are unknown. Results 
from several experiments strongly suggest that neutrino 
flavor mixing occurs in solar, atmospheric, and acceler- 
ator neutrinos, and proof of mixing would be a proof of 
mass. Direct searches for neutrino mass yield only the 
following hmits: m^^ < 5 eV [|, m,^ < 170 keV §, and 
m^^ < 24 MeV ||]. With current techniques, it will be 
very difHcult to significantly improve these limits. In fact, 
the interesting mass scale is much lower than the latter 
two limits, and is given by the requirement that neutri- 
nos do not overdose the universe (see and references 
therein): 



3 



4=1 



< 100 eV . 



(1) 



Neutrino masses exceeding this bound are allowed for 
unstable neutrinos. 

The most promising technique for direct determination 
of neutrino mass below the cosmological bound seems 
to be from time-of-flight measurements over astrophysi- 
cal distances. With the present generation of detectors, 
neutrinos and antineutrinos of all flavors from a Galactic 
supernova will be readily detectable. Even a tiny mass 
will make the velocity slightly less than for a massless 



neutrino, and over the large distance to a supernova will 
cause a measurable delay in the arrival time. A neutrino 
with a mass m (in eV) and energy E (in MeV) will ex- 
perience an energy-dependent delay (in s) relative to a 
massless neutrino in traveling over a distance D (in 10 
kpc) of 



At(£;) = 0.515 (-) D, 



(2) 



where only the lowest order in the small mass has been 
kept. For a v^^ or Vr mass near the cosmological bound, 
the delay for a single neutrino is of order seconds. Be- 
cause of the limit on the De mass, this delay can be mea- 
sured from the arrival of the and events relative 
to the Df, events. With the statistical power of many 
events, it is possible to detect an average delay of order 
0.1 s. Since one expects a type-II supernova about every 
30 years in our Galaxy , and since supernova neutrino 
detectors are currently operating, there is a good chance 
that this technique can be used to dramatically improve 
the limits on the and Vr masses. 

In a previous paper 1^, we considered a future su- 
pernova at 10 kpc (approximately the distance to the 
Galactic center) as seen by the SuperKamiokande (SK) 
detector. Such a supernova will cause about 710 neutral- 
current excitations of ^^O by and Vr (and their an- 
tiparticles), followed by detectable gamma emission. In 
addition, about 8300 events are expected from j/g -f p — > 
e^+n. A Vfj^ or Vr mass would cause a delay of the average 
arrival time of the neutral-current events as compared to 
the Ve events. We have shown how to test the statistical 
signiflcance of the difference in average arrival times and 
how to extract the allowed neutrino mass range. Taking 
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into account the finite statistics, we concluded that with 
this signal at SK, one can reach a mass sensitivity down 
to about 45 eV for the or Vt mass. 

In this paper, we consider the capabilities of the Sud- 
bury Neutrino Observatory (SNO). For the same super- 
nova, SNO will see in total about 400 events caused by 
the neutral-current breakup of deuterons by v^i and Vr 
(and their antiparticles). As in Ref. 1^, the technique 
used is to compare the average arrival time of the (pos- 
sibly massive) and Vt events with the average arrival 
time of the events. The sensitivity of SNO for this 
measurement has been estimated previously with 
the claimed minimal detectable mass ranging from 10 
eV to 200 eV. Here, we present a detailed calculation 
of the mass sensitivity of SNO, taking the finite statis- 
tics into account quantitatively. While the statistics are 
lower than for SK, the characteristic energy is lower (and 
so the delay is larger) , leading to a sensitivity to a or 
j^T- mass down to about 30 eV. 

The problem of or mass determination with su- 
pernova neutrinos has been discussed for other neutrino 
detection reactions and using different analysis tech- 
niques. Neutrino-electron scattering in water-Cerenkov 
detectors (e.g., SK) has been considered in Refs. [p|-p^. 
There have also been proposals to use neutral-current ex- 
citation of various nuclei as a signal ||l^-|l7|l. As shown 
below, some of the previous estimates seem to be too 
optimistic in that they assume a very sharp pulse of neu- 
trinos in time, which makes any delay more apparent. 

In Section II, we describe the details of the supernova 
model, the detector properties, and the neutrino detec- 
tion signals. In Section III, we review our analysis tech- 
nique and show our results for the sensitivity of SNO 
to small or Vr masses. We also consider how the 
mass sensitivity would be modified if the actual super- 
nova parameters differ from those assumed here. In Sec- 
tion IV, we discuss how the supernova parameters and 
neutrino properties can be separately extracted from the 
same data. In Section V, we summarize our results. 

II. PRODUCTION AND DETECTION OF 
SUPERNOVA NEUTRINOS 

A. Supernova neutrinos 

When the core of a large star (Af > SAfo) runs out 
of nuclear fuel, it collapses and forms a proto-neutron 
star with a central density well above the normal nu- 
clear density (for a review of type-II supernova theory, 
see Ref. Q). The total energy released in the col- 
lapse, i.e., the gravitational binding energy of the core 
(Eb ~ GnMI/R with i? - 10 km), is about 3 x 10^^ 
ergs; about 99% of that is carried away by neutrinos and 
antineutrinos, the particles with the longest mean free 
path. The proto-neutron star is dense enough that neu- 
trinos diffuse outward over a timescale of several seconds. 



maintaining thermal equilibrium with the matter. When 
they are within about one mean free path of the edge, 
they escape freely, with a thermal spectrum characteris- 
tic of the surface of last scattering. The luminosities of 
the different neutrino flavors are approximately equal at 
all times; see e.g., Ref Q. 

Those flavors which interact the most strongly with the 
matter will decouple at the largest radius and thus the 
lowest temperature. As explained in Ref. the and 
Vr neutrinos and their antiparticles, which we collectively 
call Vx neutrinos, have a temperature of about 8 MeV (or 
{E) ~ 25 MeV). The neutrinos have a temperature of 
about 5 MeV {{E) ~ 16 MeV), and the neutrinos 
have a temperature of about 3.5 MeV {{E) ~ 11 MeV); 
see e.g., Refs. fisl , ^ . These are the temperatures used 
in our analysis. While there is some variation between 
models in the actual values of the temperatures, all of 
the models have a temperature hierarchy as above. This 
is important for separating the Ux neutrinos from the 
and Df, neutrinos. The energy distributions are taken 
here to be Fermi-Dirac distributions, characterized only 
by the temperatures given above. More elaborate models 
also introduce a chemical potential parameter to reduce 
the high-energy tail of the Fermi-Dirac distribution; the 
effect of this is considered below. 

While some numerical supernova models have temper- 
atures decreasing with time, more recent models p9[ | have 
temperatures increasing with time. This is a consequence 
of the electron fraction and hence the opacities decreas- 
ing with time. The real temperature variation is probably 
not large (see e.g., Ref. ||l^]). A well- motivated form for 
temperature variation may eventually be obtained from 
the supernova D^, data or from more-developed numeri- 
cal models. The analysis of this paper could be easily 
modified to allow a varying temperature; until there is 
a compelling reason to use a particular form, we simply 
use constant temperatures. 

The neutrino luminosity rises quickly over a time of 
order 0.1 s, and then falls over a time of order several 
seconds. The luminosity used here is composed of two 
pieces. The first gives a very short rise from zero to 
the full height over a time 0.09 s, using one side of a 
Gaussian with a = 0.03 s. The rise is so fast that the 
details of its shape are irrelevant. The second piece is 
an exponential decay with time constant r = 3 s. The 
luminosity then has a width of 10 s or so, consistent with 
the SN 1987A observations. The detailed form of the 
neutrino luminosity is less important than the general 
shape features and their characteristic durations. 

This description of a supernova is consistent with the- 
oretical expectations, numerical supernova models, and 
the SN 1987A observations. With the next Galactic su- 
pernova, there will obviously be great improvements in 
the understanding of the supernova neutrinos. In Section 
IV, we discuss how to separately extract the supernova 
parameters and neutrino properties from the same data. 
Throughout the paper, we assume that the distance to 
the supernova is D = 10 kpc. 
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B. General form of the neutrino scattering rate 

Here we briefly summarize the notation of Ref. 
Under the assumption that the neutrino energy spec- 
tra are time-independent, the double differential num- 
ber distribution of neutrinos of a given flavor (one of 
j/g, i/g, t'/j, f'/j, I'r , j^r) at the source can be written as: 



dEdU 



f{E) 



L{U) 
(E) 



(3) 



where E is the neutrino energy and ti is the emission 
time. Here f{E) is the normalized thermal spectrum, 
L{ti) is the the luminosity (energy flux per unit time), 
and (E) is the (time-independent) average neutrino en- 
ergy. The double integral of this quantity is the total 
number N^, of emitted neutrinos of that flavor. This form 
is convenient since we assume that the luminosities of the 
different flavors are approximately equal at every time ti. 
As stated above, the energy spectrum f{E) is assumed to 
be a Fermi-Dirac distribution, and the luminosity L{ti) 
is assumed to have a very sharp rise and an exponential 
decline. 

The arrival time of a neutrino of mass m at the detector 
is, t = ti + D + At{E), where D is the distance to the 
source, and the energy-dependent time delay is given by 
Eq. (||). For convenience, we drop the constant D. Then 
the double differential flux of neutrinos at the detector is 
given by 



1 d^N^ 
AnD^ dEdt 



1 



1 



dt, 



d^N, 



f ,5(t - - At{E)) 



fiE) 



dEdt, 
L{t~ At{E)) 



(4) 



Note that because of the mass effects, this is no longer 
the product of a function of energy alone and a function 
of time alone. The scattering rate for a given neutrino 
reaction in SNO is 



dN^ 
dt 



Nb^o n / dE(j{E) 



1 d^N^ 
47rL>2 dEdt 



(5) 



where Nd^^q is the number of heavy- water molecules in 
the detector, a{E) the cross section for a neutrino of 
energy E on the target particle, and n the number of 
targets per molecule for the given reaction. Using the 
results above, the scattering rate (in s~^) can be written: 



dN^ 
dt 



C = 8.28 



C / dEf{E) 



a{E) 



10- 



L{t-At{E)) 



Eb/6 



(6) 



Eb 



IQ53 gj.gg 



1 MeV 



T 



10 kpc 



D 



det. mass 
1 kton 



(7) 



In the above, T is the spectrum temperature (where we 
assume (E) = 3.15 T, as appropriate for a Fermi-Dirac 
spectrum), and f{E) is in MeV^^. For a light- water 
detector, the initial coefficient in C is 9.21 instead of 8.28. 
For equal luminosities in each flavor, the total binding 
energy released in a given flavor is Eb/6- Note that we 
ignore the prompt burst of Ve neutrinos, since these carry 
only of order 1% of the total energy. When an integral 
over all arrival times is made, the luminosity term in 
Eq. (^) integrates to one, giving for the total number of 
scattering events: 



N,. 



C / dEf{E) 



a{E) 



10- 



(8) 



For massless neutrinos, At{E) — 0, and in Eq. (||) the 
luminosity can be taken outside of the integral, so that 
the time dependence of the scattering rate depends only 
on the time dependence of the luminosity. 



C. The Sudbury Neutrino Observatory 

When it begins operations this year, the Sudbury Neu- 
trino Observatory (SNO) will be the first deuterium- 
based detector for neutrino astrophysics. The SNO de- 
tector is described in Refs. [0|2^. Here we give a short 
summary of the properties relevant to our analysis. Deu- 
terium is an excellent target for neutrinos since both the 
charged- and neutral-current cross sections are reason- 
ably large. The active part of the detector is 1 kton of 
pure D2O, separated by an acrylic vessel from 1.7 kton 
of light water. This entire volume is viewed by 10** pho- 
totubes which can see about 1.4 kton of the light water 
with good efficiency We assume that events in the 
D2O can be separated from events in the H2O. 

The proposed threshold for solar neutrino studies is 5 
MeV (the physics potential of the SNO solar neutrino 
studies is treated in Ref. |^ and references therein). At 
this energy, the contribution from the time-independent 
background is expected to be small [|l|. For a Galactic 
supernova, one expects several hundred events over about 
10 s, and a much higher background rate can be tolerated, 
allowing a lower threshold. For a threshold of 5 MeV, the 
solar neutrino rate of about lO^** s"^ is a background for 
supernova neutrinos. Below 5 MeV, the background rate 
increases very steeply. From Ref. |2^, we estimate that 
the threshold for the supernova analysis can be lowered 
by a few MeV while keeping the background contribution 
negligible. This would ensure that almost all low-energy 
gammas from neutron captures as well as electrons and 
positrons from charged-current events are detected. 

Electrons and positrons will be detected by their 
Cerenkov radiation, and gammas via secondary electrons 
and positrons. It is not possible for SNO to distinguish 
between electron, positrons, and gammas of compara- 
ble energy. The only way to detect the neutral-current 
breakup of the deuteron is to detect the final neutron. 
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There are three neutron detection techniques proposed 
for SNO: (n,7) on deuterons in pure D2O, giving a 
gamma of energy 6.25 MeV; {n,j) on ^^Cl in D2O with 
MgCl2 sah added, giving a gamma cascade of energy 8.6 
MeV; and direct n detection in ^He proportional counter 
tubes suspended into the D2O. If either of the latter two 
techniques are used, there will still be some contribution 
from neutron captures on deuterons. In this paper, as in 
all previous studies of the supernova capabilities of SNO, 
we assume that the neutron detection efficiency is nearly 
100%. In Section III, we show that a reduced neutron 
detection efficiency has only a small effect on the mass 
sensitivity. It therefore makes very little difference which 
neutron detection technique or techniques arc in place 
when the supernova occurs. 

D. Neutrino signals in SNO 

Neutrinos and antineutrinos of all flavors cause the 
neutral-current breakup of the dcutcron: 

v + d^v+p + n, (9) 

D + d^ D + p + n, (10) 

with thresholds of 2.22 MeV. There are about 35 Ve, 110 
j/^, and 110 i>r events from the first reaction, and 50 
j>e, 90 Dfi, and 90 Dr events from the second reaction. 
Electron neutrinos and antineutrinos in addition cause 
the charged-current breakup of the deuteron: 

i^e + d~^e^+p + p, (11) 

i?e + d^e+ + n + n, (12) 

with thresholds of 1.44 and 4.03 MeV, respectively. There 
are about 80 events each from these two reactions. The 
numbers of events are summarized in Table I. The cross 
sections for all these reactions have an energy dependence 
roughly of the form cr{E) ^ {E — £'th)^- There have 
been many calculations of the neutrino-deuteron cross 
sections |p^ , p5| , and they are now rather well-determined. 
In this paper, we use the tabulated cross sections of 
Ref. Il^; the quoted uncertainty is 5% or less at the 
relevant energies. All of the outgoing particles in the 
deuteron breakup reactions are emitted approximately 
isotropically. 

The neutral-current reactions are flavor-blind. How- 
ever, since neutrinos and antineutrinos from a supernova 
have spectra with a hierarchy of temperatures, the energy 
dependence of the rates favors the higher-temperature 
flavors. In particular, most (about 82%) of the events 
will be from i^^ neutrinos. (Previous studies which indi- 
cated a percentage of about 90% used a higher tem- 
perature.) 

In order to estimate the delay, Eq. (|^) can be evaluated 
with a characteristic energy. However, one should not 



TABLE I. Calculated numbers of events expected in SNO 
for a supernova at 10 kpc. The other parameters (e.g., neu- 
trino spectrum temperatures) are given in the text. In rows 
with two reactions listed, the number of events is the total 
for both. The notation 1/ indicates the sum of Ue, v^, and 
Vt, though they do not contribute equally to a given reaction, 
and X indicates either or p-(-^^N. 



Events in 1 kton D2O 



V -\- d ^ ly + p + n 


485 


V + d ^ iy + p + n 




Ue + d^e^+p + p 


160 


Ue + d—>e^+n + n 




iy+"'0^v + 'y + X 


20 


u+^'^O^u + 'y + X 






15 






u + e~ 1/ + 


10 


V + e~ + 




Events in 1.4 kton H2O 


Ue + p ^ e''' + n 


365 


u+'-'^O^u + 'y + X 


30 


v+^'^O^u + j + X 




1/ + e~ f + 


15 


u + u + 





use the average energy (E), but rather a characteristic 
energy that takes into account the weighting with the 
cross section as well. For the neutral-current neutrino- 
deuteron reactions, this is shown later to be about 32 
MeV. In addition, the fact that the neutrinos have a 
spectrum of energies means that different values of E 
contribute to the time delay, causing dispersion of the 
neutrino pulse as it travels from the supernova. It turns 
out that for the small masses we are primarily interested 
in, these dispersive effects are minimal. 

Since the neutron capture time may be as large as sev- 
eral ms, and the event rates may be high, concerns were 
raised in Refs. [|lO|,0 that events would overlap in time 
and that it would therefore be difficult to distinguish 
charged-current and neutral-current events. Only dur- 
ing the first second or so are the rates likely to be high 
enough that this can occur. In the first second, there 
are about 150 neutral-current events and only about 40 
charged-current events. Further, if the neutron mean free 
path is less than the diameter of the acrylic vessel, then 
it will be possible to use spatial information to distin- 
guish events. In any case, the possible contamination 
of the neutral-current rate is very small, and the mass 
sensitivity will not be significantly affected. 

For the solar neutrino studies, the electron from the 
charged-current reaction has a low energy and can be 
confused with a gamma from a neutron-capture event. 
Since the supernova i^e energy is much higher, the elec- 
tron energy in a charged-current reaction is high enough 
that only rarely can it be confused with a gamma. Thus 
neutrons can almost always be identified, either by direct 
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detection with a proportional counter, or by the energy of 
the subsequent gamma. Positrons from charged-current 
events with Ve can be identified by their high energy and 
coincidence with two neutrons. The spectra of electron 
and positron energies will be broad, peaking at about 15 
and 20 MeV, respectively. 

So far, we have discussed only the neutrino-deuteron 
reactions. As noted in Ref. |^^, the neutral-current ex- 
citation of ^^O into the continuum, followed by neutron 
or proton emission, can also be detected. About 30% of 
the time, the A = 15 nucleus is left in an excited state 
which decays by gamma emission, with gamma energies 
between 5 and 10 MeV. These gammas are detectable, 
as is the neutron. The remaining 70% of the time, the A 
— 15 nucleus is left in the ground state. In these cases, 
only the final states with a neutron are detectable. All of 
the outgoing particles are emitted approximately isotrop- 
ically. Including the "'^^O events only has a small effect 
on our final results. 

There are also events from neutrino-electron scatter- 
ing, which we ignore. These are forward-peaked, and we 
assume that they have been removed with an angular 
cut. A forward cone of half-angle 25 degrees would con- 
tain almost all neutrino-electron scattering events, while 
removing only 5% of the isotropic events. 

III. SIGNATURE OF A SMALL NEUTRINO 
MASS 

A. General description of the data 

For a massless neutrino {v^. or v^) the time dependence 
of the scattering rate is simply the time dependence of the 
luminosity of that flavor. For a possibly massive neutrino 
{vx), the time dependence of the scattering rate depends 
not only on the time dependence of the luminosity of that 
flavor, but also on the delaying effects of a mass. Since 
the luminosities of the different flavors are expected to 
be equal at all times, then if the is massive, we can 
compare the Vx and Ve scatte ring rates to search for a 
Vx mass; see Eq. d). In Rcf. Q, it was proposed that 
the effects of a mass could be determined from the time 
dependence of the neutral-current rate divided by the 
total rate, which also has the effect of removing the the 
time dependence of the luminosity from the scattering 
rate. However, no quantitative technique for extracting 
the allowed mass range was given. 

In order to implement the comparison of scattering 
rates, we define two rates: a Reference R{t) containing 
only massless events, and a Signal S{t) containing some 
fraction of possibly massive events, and some fraction 
of massless events (the possibly massive events cannot 
be completely separated from the massless events). We 
assume that only Vt is massive (this cannot be distin- 
guished from the case that only is massive). The 
analysis can easily be repeated for the case that both 



and Vt are massive. Because of the quadratic depen- 
dence of the delay on the mass, this two-mass case will 
look like a one-mass case with the larger mass unless the 
two masses are nearly equal. 

The Reference R{t) can be formed in various ways. The 
largest sample of useful massless events will be the 7800 
De events from SK with > 10 MeV. Below 10 MeV, 
there are gammas from the neutral-current excitation of 
^^O which cannot be separated. We assume that SK and 
SNO will have synchronized clocks so that in principle, 
such a sharing of data will be possible. One can also 
use the 340 Vf, events with Ef^+ > 10 MeV in the light 
water at SNO. Because of the smaller number of counts, 
the statistical error is larger and hence the mass sensi- 
tivity is slightly worse. This could be slightly improved 
by including the 160 charged-current events in the heavy 
water (a small fraction of the low-energy events would 
again have to be cut). 

The primary component of the Signal S{t) is the 485 
neutral-current events on deuterons. With the temper- 
atures assumed here, these events are 18% {vg + i^e), 
41% (f^ -I- Df^), and 41% + Vr). The flavors of the 
neutral-current events of course cannot be distinguished. 
Therefore, under our assumption that only Vr is mas- 
sive, there is already some unavoidable dilution of the ex- 
pected delaying effect of a mass. The Signal S{t) should 
also contain all events below about 10 MeV which can- 
not otherwise be removed. There are 35 neutral-current 
events on ^^O in the heavy water. Because of the high 
threshold for this reaction, these events are 50% each of 
{v^ -\- v^), and {vt + P^). Finally, low-energy charged- 
current events must also be included. Because of the 
high energies of supernova neutrinos, only a small num- 
ber (about 15) of electrons must be included in the Signal 
S{t). No positrons need be included since those events 
can be identified by their two accompanying neutrons. 
The inclusion of the ^^O and charged-current events only 
slightly changes the shape of S{t), and barely changes the 
final mass limit. The total number of events in S{t) is 
535, of which still 41% are possibly massive. 

If the Signal S{t) contains some fraction of massive 
events, the shape of the scattering rate will be delayed 
and broadened. Relative to the Reference R{t), there is 
a deficit of events at early times and an excess at late 
times. In the test for a mass, we test for this characteris- 
tic distortion in the shape. In Fig. 1, S{t) is shown under 
different assumptions about the v-r mass. The shape of 
R{t) is exactly that of S{t) when m^^ = 0, though the 
number of events in R{t) will be greater than in S{t) if 
the SK R{t) is used and comparable if the SNO R{t) is 
used. For a very large v-r mass, the massless and massive 
components of S{t) would completely separate in time. 
In Fig. 1, note that for m — 100 eV, almost all of the 
massive events are delayed beyond 1 s. 

The rates R{t) and S{t) will be measured with finite 
statistics, so it is possible for statistical fluctuations to 
obscure the effects of a mass when there is one, or to fake 
the effects when there is not. From Fig. 1, and Poisson 
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FIG. 1. The expected event rate for the Signal S{t) at SNO 
in the absence of fluctuations for different i/j- masses, as fol- 
lows: solid line, eV; dashed lines, in order of decreasing 
height: 20, 40, 60, 80, 100 eV. Of 535 total events, 100 are 
massless {i/^ + v^), 217.5 are massless (i/^ + i/p), and 217.5 are 
massive {vr + Vt)- These totals count events at all times; in 
the figure, only those with i; < 9 s are shown. 



statistics, one can easily get a rough idea of how finely 
the mass can be determined from the difference between 
R{t) and S{t). Note that if the SK R{t) is used, the 
fluctuations in R{t) when scaled down to the number of 
events in S{t) will be small, and that if the SNO R{t) is 
used, the fluctuations in R{t) will be comparable to those 
in S{t). 

In this paper, we determine the mass sensitivity in the 
presence of the statistical fluctuations by Monte Carlo 
modeling. We use the Monte Carlo to generate represen- 
tative statistical instances of the theoretical forms of R{t) 
and S{t)^ so that each run represents one supernova as 
seen in SNO. For R{t), we pick a Poisson random num- 
ber from a distribution with mean given by the expected 
number of events. This determines the number of events 
for a particular instance. We then use an acceptance- 
rejection method to sample the form of R{t) until that 
number of events is obtained. This gives a statistical in- 
stance of representative of what might be seen in a 
single experiment. A similar technique is used to gener- 
ate an instance of S{t). The massless and massive com- 
ponents of S{t) are sampled separately, and then added 
together. 

In Ref. |q|, we considered two different tests of the 



shape distortion of S{t) relative to R{t). The first was a 
test. A large value of between statistical instances 
of S{t) and R{t) would indicate that they were likely to 
have been drawn from different distributions, and this 
would be taken as evidence of a v-r mass. However, this 
test is non-specific to testing for a mass, i.e., it is sensitive 
to any difference between S{t) and R{t). In addition, the 
test requires binning in time, which washes out the 
effects of small masses. The second test, and the one ad- 
vocated there, was a test of the average arrival time (i). 
Any massive component in S{t) will always increase (i), 
up to statistical fluctuations. Besides being more directly 
related to the mass effect, the {t) technique is sensitive to 
somewhat smaller masses than the technique, since no 
binning is required. In order that the results be believ- 
able, it is necessary that different reasonable statistical 
techniques yield consistent results. For this paper, we 
used both techniques and verified that they gave similar 
results. However, we present only the results of the (t) 
analysis. 



B. {t) analysis 

Given the Reference the average arrival time is 

defined as 



{t)i 



(13) 



The summation form is used for real or simulated data 
sets, where the sum is over events (not time bins) in the 
Reference with < t < t^ax- The integral form would 
be used if the theoretical form for the rate were given. 
The starting time is assumed to be well-defined. With 
some 10^ events expected in SK and SNO combined, and 
a risetime of order 0.1 s, this should not be a problem; 
the definition used here amounts to calling the starting 
point that time at which the J>e rate is about 1% of its 
peak rate. The choice of tmax is made as follows. The 
effect of the finite number of counts in R{t) is to give {t)ii 
a statistical error: 



(14) 



where both the width {t"^) n — (t)^ and the number of 
events Nr depend on t,nax- By choosing a moderate 
tmaxi the width of R{t) can be restricted. Such a choice 
will make the error on as small as possible given 
that almost all of the events are to be included. Both 
the starting time and tmax were of course held constant 
over different Monte Carlo runs. For a purely exponential 



luminosity, and tr, 



oo, {t)R = 



Given the Signal S{t)^ the average arrival time is de- 
fined similarly as 



(15) 



6 



where naturally the sums are now over events in the Sig- 
nal. The widths of R{t) and S{t) are similar, each of 
order r = 3 s (the mass increases the width of S{t) only 
slightly for small masses.) If the SK R{t) is used, then 
the statistical error on (t) s is few times larger than that 
on {t)ii since there arc several times fewer events. If the 
SNO R{t) is used, then the statistical errors on {t)s and 
{t)R are comparable. Note that the errors on {t)^ and 
{t)s are uncorrelated. 

The signal of a mass is that the measured value of 
{t)s — {t)R is greater than zero with statistical signifi- 
cance. From the Monte Carlo studies, tmax = 9 s was 
found to be a very reasonable choice for the luminos- 
ity decay time r = 3 s; about 95% of the data are 
then included while the width is somewhat reduced. The 
time-independent background event s are negligi ble. For 
tmax = 9 s, {t)R = 2.57 s and ^{t^)R - {t)% = 2.12 
s. Near tmax = 9 s, the significance of a delay, i.e., 
{t)s ~ {t)R divided by its statistical error, is nearly max- 
imal for the small masses wc are considering. However, 
the results are not strongly dependent on the particular 
value of tjnax used as long as it is reasonable (i.e., the vast 
majority of events are contained) . Note that any shift in 
the starting time will cancel in the difference {t)s — {t)R- 

Using the above procedure, we analyzed 10^ simulated 
supernova data sets for a range of Vr masses. For each 
of them, {t)s — {t)R was calculated and its value his- 
togrammed. These histograms are shown in the upper 
panel of Fig. 2 for a few representative masses. (Note 
that the number of Monte Carlo runs only affects how 
smoothly these histograms arc filled out, and not their 
width or placement.) These distributions are character- 
ized by their central point and their width, using the 
10%, 50% (equal to the average), and 90% confidence 
levels. That is, for each mass we determined the values 
of {t)s — {t)R such that a given pcax-cntagc of the Monte 
Carlo runs yielded a value of {t)s — {t)R less than that 
value. With these three numbers, we can characterize the 
results of complete runs with many masses much more 
compactly, as shown in the lower panel of Fig. 2. Since 
the {t)s — distributions are Gaussians, other confi- 
dence levels can easily be constructed. For convenience, 
the axes in the lower panel are inverted from how the 
plot was actually constructed. That is, given an experi- 
mentally determined value of (t) s — (t) r, one can read off 
the range of masses that would have been likely (at these 
confidence levels) to have given such a value of {t)s — {t)R 
in one experiment. From the lower panel of Fig. 2, we see 
that SNO has a sensitivity to a Vr mass down to about 
30 eV (rounded from 27.5 cV) if the SK R{t) is used, and 
down to about 35 eV if the SNO R{t) is used. 

We also investigated the dispersion of the event rate 
in time as a measure of the mass. A mass alone causes 
a delay, but a mass and an energy spectrum also cause 
dispersi on. We define d t he dispersion as the change in the 
width [t"^) s — — \/ (i^)fl — where all integrals 
are as above defined up to tmax- We found that the 
dispersion was not statistically significant until the mass 
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<t>s-<t>^ [s] 
FIG. 2. The results of the {t) analysis for a massive v-r, 
using the Signal S(t) from SNO defined in the text. In the 
upper panel, the relative frequencies of various (t)s — {t)R val- 
ues are shown for a few example masses. The solid line is for 
the results using the SK Reference R{t), and the dotted line 
for the results using the SNO R{t). In the lower panel, the 
range of masses corresponding to a given {t)s — {t}R is shown. 
The dashed line is the 50% confidence level. The upper and 
lower solid lines are the 10% and 90% confidence levels, re- 
spectively, for the results with the SK R{t). The dotted lines 
are the same for the results with the SNO R{t). In this fig- 
ure, tmax = 9 s and the time constant of the exponential 
luminosity is t = 3 s. 



was of order 80 eV or so; however, for such a large mass 
the statistical significance of the change in (t) cannot be 
missed. For these large masses, dispersion docs increase 
the width of S{t) and hence the error on {t)s — {t)R; this 
is just becoming visible in Fig. 2. 



C. Sensitivity to the input parameters 

Since the parameters governing the supernova neutrino 
emission are not perfectly known, it is worthwhile to 
examine the sensitivity of our conclusions to their as- 
sumed values. Most of this dependence can be obtained 
analytically by extending the time integration limit to 
tmax oo. Note that the mass sensitivity will be slightly 
poorer than for tmax = 9 s used in the main analysis. 

The characteristic delay is then 
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{t)s - {t)R ^ frac(TO > 0) X 0.515 { — \ D. 



(16) 



where frac(m > 0) is the fraction (about 41%) of mas- 
sive events in S(t) and the units are as in Eq. (H). The 
characteristic energy Ec can be taken to be the peak of 
f{E)a{E), or more accurately from 



Er. 



JdEf{E)a{E){^ 
JdEfiE)a{E) 



(17) 



The delay corresponding to Ec is the delay of the centroid 
of S{t) relative to the centroid of R{t), i.e., exactly {t)s — 
{t)n- For SNO, Ec ~ 32 MeV. The statistical error on 
{t)ii is given by Eq. (^4|), and the statistical error on (t)s 
is defined similarly. 

Using the above, we can make a good estimate of the 
mass limit that SNO would set if the delay were measured 
to be zero. Ignoring the error on (i)/?, and taking t^ax ~^ 
oo, we estimate the error to be t/^Ns — 3/\/535 — 0.13 
s. The 10% and 90% confidence levels as used in Fig. 2 
correspond to ±1.3 x 5 {{t)R — {t)s)- For our estimate, 
this has magnitude 0.17 s, close to the result in Fig. 2 for 
the SK Reference. The mass limit that will be placed if 
no delay is seen (using D = \, i.e., 10 kpc) is 



'^livait — -^c 



1.3 X error 



frac(m > 0) x 0.515D 



30 eV, (18) 



in excellent agreement with the numerical result. If the 
neutron detection efficiency e„ were not 100%, then Ng 
would be reduced by Cn and the error increased by 1/ 
Therefore, the mass limit would be increased only by 

This formula can also be used to make an estimate of 
the mass limit that can be obtained with the forward- 
peaked neutrino-electron scattering at SK (which has 
the largest number of those events). Since (j{E) ^ E, 
Ec ~ {E) ~ 25 MeV. The Signal S{t) can be defined by 
the events in a forward cone of half-angle 25 degrees, 
and the Reference R{t) defined by the events outside 
this cone. Assuming that all of the neutrino-electron 
scattering events are contained along with about 5% 
of the isotropic backgrounds, Ng = 760 and frac(m > 
0) = 60/760 = 0.079 [|. The error is estimated to be 
3/v'760 — 0.11 s, and so mumit — 50 eV. This is smaller 
than the recent conclusion of Ref . , but the technique 
used here includes a much greater portion of the data. 

Continuing to suppose that the mass is very small, and 
that SNO will simply place a limit, we can investigate 
the effects of varying the input parameters. We ignore 
dispersion and take the widths of R{t) and S{t) to be 
proportional to r. If the cross section (t{E) depends on 
energy as i?" (a ~ 2 for + d), then the characteristic 
energy Ec ^ [2 + a)T and the thermally-averaged cross 
section is proportional to T", where T is the tempera- 
ture. We take frac(m > 0) to be approximately constant 



and ignore the error on (t)_R; this is only valid for small 
deviations from T = 8 MeV. With these assumptions we 
can determine how the key parameters affect the result 
(for a = 2). The delay is 



{t)s - {t)R 

The number of events is 



/m 



D. 



T 



and so the error is 

H{t)s-{t)R) 



(19) 



(20) 



(21) 



Therefore, the significance (number of sigmas) of the de- 
lay is 



{t)s - {t)R 

sms~{t)R) 



T 



T 



(22) 



Remarkably, this is independent oi D. To place a mass 
limit at a given confidence level, the number of sigmas 
is fixed. Then the mass limit that can be obtained if no 
delay is observed has the following dependence: 



rp3/4 



(23) 



also independent of D. 

The D-independence can also be seen geometrically 
from Fig. 2. Under D 2D, for example, the delay and 
the error each become twice as large. However, the point 
at which the upper confidence level crosses the m-axis 
is unchanged. If a nonzero delay is measured, then the 
range of allowed masses does depend on D, with larger 
distances corresponding to smaller masses. The above ar- 
guments hold for the relevant range 1-30 kpc; for smaller 
D the detector may be overwhelmed and for larger D 
there are too few events and the derivation is not valid. 
Because of obscuration by dust, it may be difficult to 
measure the distance to a future Galactic supernova. It 
is therefore rather important that this does not affect 
the ability to place a limit on the Vr mass. That may 
not be true for other analysis techniques; for example, 
the estimates of Ref. Q do depend on D. Note that a 
smaller distance would allow a better measurement of the 
temperatures and other supernova parameters. 

Now consider the effect of a change in the temperature 
on the mass limit. As indicated, the dependence is 
weak, with a higher temperature making it slightly more 
difficult to limit the mass. The numerical results and the 
analytic estimate are shown in Fig. 3. Even under the 
large variation in T,y^ of ±2 MeV, the mass limit changes 
only by about ±5 eV. 

The effect of changing the luminosity decay time con- 
stant T is even more straightforward and is obvious from 
Eq. (p3[). Thus, reducing that from T = 3stoT = l 
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FIG. 3. The effect of a change in the temperature T 
on the mass limit that could be made using the SNO Sig- 
nal and the SK Reference. The mass limit is defined as the 
mass at which the appropriate upper confidence level inter- 
sects the m-axis in Fig. 2. The solid line is the formula given 
in Eq. (^). The full numerical results are marked by circles, 
and connected with a dashed line to guide the eye. 

s, as used in Ref. [|o|, would improve the limit at SNO 
(using the SK Ref) from about 30 eV to about 15 eV. 
Using T — 0.5 s, as in Ref. ||l^, or r = 0.3 s, as in 
Ref. nj, would improve the limit to about 10 eV. These 
values of t were estimated from appropriate figures in 
these references, and are only approximations of the L(t) 
timescale. The numerical results and the analytic esti- 
mate are shown in Fig. 4. Using an unrealistically sharp 
neutrino pulse makes the quoted mass limit very small. 

We also considered the effect of a chemical potential 
in the spectrum, which would reduce the high-energy 
tail. We took ^ = 3T, and then chose T = 6.31 MeV to 
keep the average energy the same as for the /i = 0, T = 8 
MeV case. Because the neutrino-deuteron cross section 
only depends quadratically on energy, the effect is small. 
About 10% fewer events are obtained, but with a delay 
about 10% larger. These lead to a change of order a few 
percent in the mass sensitivity. Because of their steeper 
energy dependence the i^x+^^O events would be more 
affected by a chemical potential. 



I • ■ • — • • ■ — I 

0.1 1.0 

T [S] 

FIG. 4. The effect of a change in the luminosity decay con- 
stant r on the mass limit that could be made using the SNO 
Signal and the SK Reference. The mass limit is defined as 
the mass at which the appropriate upper confidence level in- 
tersects the m-axis in Fig. 2. The solid line is the formula 
given in Eq. (^3|). The full numerical results are marked by 
circles, at the points r = 3 s (this paper), r = 1 s (Ref. ]lo|]), 
T = 0.5 s (Ref. Q), and r = 0.3 s (Ref. 0); see the text 
for explanation. 

IV. SEPARATE EXTRACTION OF SUPERNOVA 
AND NEUTRINO PARAMETERS 

The observation of the neutrino signal of a future 
Galactic supernova will be extremely significant test of 
the physics involved. It will allow, among other things, 
determination of the imprecisely-known supernova neu- 
trino emission parameters. In addition, we hope to be 
able to use the same data to determine or constrain neu- 
trino properties. In this section, we discuss how both of 
these goals can be achieved simultaneously. Throughout, 
we use standard values for all parameters to numerically 
evaluate the expected precision. 

The measured neutrino signal can be used to deter- 
mine the temperatures (and possibly the chemical po- 
tentials) of Ve, f'e, and Vx- Ideally, we would like to 
know these parameters as a function of time. In addi- 
tion, the data can be used to determine the apparent 
source strength Eb /D^ and the luminosity time profile. 
Since neutrino oscillations could change the neutrino fla- 
vor between emission and detection, their presence would 
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make this task more difficult. Thus, we consider first the 
simpler situation without neutrino oscillations. 

1. Without neutrino oscillations 

The Vf, and neutrinos will be detected by the 
charged-current reactions in which the energy of the out- 
going electron or positron will be measured. Neglect- 
ing recoil corrections and detector resolution, the relation 
El, = Ee + Eth allows one to relate the measured electron 
or positron spectrum to the weighted neutrino spectrum 
f{E^)a{E^) which appears in Eq. (||). Since the charged- 
current cross sections are well-known, the actual neutrino 
spectrum f{Eu) can be obtained. If this has the expected 
thermal shape, then the corresponding temperatures T^^ 
and Tp^ (and possibly chemical potentials) can be ex- 
tracted. The recoil corrections and detector resolution 
can of course be taken into account by fitting the mea- 
sured electron or positron spectrum with an appropri- 
ately convolved trial neutrino spectrum. Note that these 
temperature determinations depend only on the shape 
of the electron or positron spectrum. The normalization 
(the total number of events) can be used to determine 
the apparent source strength Eb/D"^. 

There are 7800 De+P ^ e+ -f n events with Ef,+ > 10 
MeV expected in SK, so that the average positron energy 
can be determined to about 1%. Following the procedure 
above, one should be able to determine the temperature 
Tp^ to a comparable precision. In this case one can also 
check for the presence of a chemical potential, and for 
time variation of Tp_, . The situation with i^e is more dif- 
ficult, since there are only 80 events expected for the 
reaction i/g + d ^ + p + p in SNO. (Note that i^e has 
no charged-current interaction with protons, that SNO is 
much smaller than SK, and that T^^ < Tp_,.) In this case, 
the average electron energy can be determined to about 
10%; as above, therefore a similar precision for the tem- 
perature T^^ . These two measurements thus allow one 
to test T^^ < Tp^ , as expected for the supernova temper- 
ature hierarchy. In addition, the two independent mea- 
surements of Eb/D^ can be tested for consistency with 
each other and with possible determinations unrelated to 
neutrinos. If Eb/D^ is assumed to be known, then T^,^ 
and Tp^ can also be determined from the total numbers 
of charged-current events in SNO, as in Ref. p7[ |. 

These measurements can also be compared with other 
available but lower-statistics data. For Tp^, other 
hydrogen-containing detectors (MACRO, SNO, Borex- 
ino, and Kamland) will have events from +p — > e+ +n. 
For T^^ , there are small numbers of events expected from 
charged-current reactions on ^^C in Borexino and Kam- 
land, and ^^O in SK. These consistency checks will be 
useful because the different reactions may have different 
systematic errors. 

The Vx neutrino can be detected only by the neutral 
current reactions. Since the outgoing neutrino carries an 



unknown amount of energy, it is not possible to deter- 
mine the neutrino spectrum as above. Therefore, T,^^ 
can only be extracted from the observed total number of 
events, and must rely on the measured value of Eb/D'^ 
and the assumption of luminosity equipartition. For the 
neutral-current reactions, the rate is independent of the 
neutrino flavor. If the supernova temperature hierarchy 
holds, then the neutral-current events will be dominated 
by the Vx neutrinos. However, one cannot directly asso- 
ciate a flavor with the extracted T^^ . 

The Signal at SK is expected to contain 710 events 
from Vx+^^O. Low-energy charged current events from 
i^e +P —^e'^+n can be confused with the neutral-current 
events; these add 530 events to the Signal. The latter 
events can be subtracted from the total number of events 
observed in the Signal since their expected number can 
be calculated using the parameters measured from the 
higher-energy De data. For the SK Signal, since the de- 
pendence of the cross section on energy is stronger, the 
4% precision on the number of events would translate to 
a 1% precision on T^,^. However, uncertainties on the 
cross section itself should also be considered. In princi- 
ple, the neutrino-electron scattering events at SK could 
also be used to measure T^^; however, the temperature 
dependence of the number of events is extremely weak. 
The Signal at SNO is expected to contain 400 events from 
Vx + d, and 135 events from the sum of neutral-current 
I'e + i'e reactions on deuterons, neutral-current excitation 
of ^^O, and low-energy charged-current reactions on 
deuterons. Again, the numbers expected for the latter 
can be calculated and subtracted from the total number 
of events observed in the Signal. The thermally-averaged 
cross section for Vx + d depends roughly on T^^ , so by 
Eq. (H) the number of events is roughly proportional to 
T^^ . Ignoring all other uncertainties, this would therefore 
allow a determination of T^^ to about 5%. There may 
also be T,^^ measurements from Vx+^^C at Borexino and 
Kamland, particularly with the excitation of the 15.11 
MeV state in ^^C. If each neutral-current reaction were 
thoroughly understood, it would be possible to use their 
different thresholds and energy dependences to map out 
the i^x spectrum. At the very least, it should be possible 
to extract a good value of T^^ and to make important 
consistency tests. 

Therefore, all three temperatures can be extracted 
with reasonable precision and the supernova tempera- 
ture hierarchy T^^ < Tp_, < T^^ experimentally verified. 
As for the time variation of the temperatures, it is likely 
that only the data will have sufficient statistics. Those 
data, after correcting for any temperature variation, will 
give the time profile of the luminosity L{t) and the char- 
acteristic timescale r. 

Mass effects will not influence the temperature deter- 
mination, provided the mass is not too large, since a de- 
lay changes only the time structure of the rate S(t), not 
its normalization nor the time-integrated energy spec- 
tra. However, if m^^ ^ 1 keV or larger, some low-energy 
events will be lost in the time-independent background 
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due to their large delay, and the extracted temperatures 
would be affected. The considerations in this section, 
combined with the result above that the mass limit is 
only weakly dependent on T^^ , shows that the temper- 
atures and the mass can be separately and robustly 
determined from the data. 



2,. With neutrino oscillations 

Assuming that the presence of oscillations does not 
substantially change the supernova dynamics, it is rela- 
tively straightforward to separately search for oscillations 
and a mass delay. The effects of supernova neutrino oscil- 
lations (without a mass delay) were considered for SNO 
in Ref. pq |. Since the numbers of i/^, j?^, i/,-, and Dr are 
expected to be the same, mixing among them have no 
effect on the number of events. What propagates are 
the mass eigenstates; the flavor content of the heavy mass 
eigenstate at the detector is irrelevant, and so oscillations 
among the neutrinos do not affect the delay. 

Now consider mixing between either j/g and Vt (or Vf^) 
or De and Dr (or P^). We assume that the Vt mass is large, 
which makes 5m^ large. Large-angle, large-^m^ mixing is 
ruled out by reactor and accelerator experiments. Small- 
angle, large-5m^ mixing is allowed; however, the effects 
are minimal unless there is an MSW enhancement. For a 
normal mass hierarchy, this can only happen for mixing 
between I'e and Vt, and at a very high density, of order 
10^^ g/cm'^, which does occur in supernovae. If there 
are such oscillations, then some high-T Vt neutrinos will 
become high-T v^. neutrinos. As noted, the spectrum of 
electrons from -\- d can be related to the spectrum of 
i^e neutrinos. For no oscillations, that would be a ther- 
mal spectrum with T^^ = 3.5 MeV; for complete oscilla- 
tions, that would be a thermal spectrum with T,^^ — 8 
MeV (and a much greater number of events). For par- 
tial mixing, there would be two peaks. The mixing pa- 
rameters and temperatures can thus be extracted from 
the measured electron spectrum. One should note that 
the charged-current reactions i^e+^^C at Borexino and 
Kamland and Ve+^^O at SK would have large numbers 
of events if such oscillations occur. If there is an in- 
verted mass hierarchy, then mixing between and 
could have an MSW enhancement. Then similar con- 
siderations to those above could be used to examine the 
positron spectrum from reactions. It would in fact be 
somewhat easier, due to the larger numbers of events. 
In either of these cases, extracting the allowed mass 
range from the measured value of (t) s — (t) r requires 
some care. A given delay could be caused by a i>r mass 
and no mixing or by a larger mass and partial mixing. 
Ideally, the mixing parameters will be known from the 
considerations above, so that the range of possible masses 
can be reasonably restricted. 



V. CONCLUSIONS AND DISCUSSION 

One of the key points of our technique is that the abun- 
dant Pe events can be used to calibrate the neutrino lumi- 
nosity of the supernova and to define a clock by which to 
measure the delay of the neutrinos. The internal cal- 
ibration substantially reduces the model dependence of 
our results, and allows us to be sensitive to rather small 
masses. Our calculations indicate that a significant delay 
can be seen for m = 30 eV with the SNO data, corre- 
sponding to a delay in the average arrival time of about 
0.15 s. Even though the duration of the pulse is expected 
to be of order 10 s, such a small average delay can be seen 
because several hundred events are expected. Without 
such a clock, one cannot determine a mass limit with the 
{t)s — {t)R technique advocated here, since the absolute 
delay is unknown. Instead, one would have to constrain 
the mass from the observed dispersion of the events; only 
for a mass of m = 150 eV or greater would the pulse be- 
come significantly broader than expected from theory. 

Moreover, the technique used here allows accurate an- 
alytic estimates of the results, so that it is easy to see 
how the conclusions would change if different input pa- 
rameters were used. If the Vr mass is very small, then the 
most probable measured delay is {t)s — {t)R = 0- In that 
case, one can only place a limit mumit on the v-r mass. 
We have shown that this limit is only weakly dependent 
on the Vx temperature T or the presence of a chemical 
potential in the thermal spectrum, and is independent 
of the distance D for the reasonable range of distances 
for a Galactic supernova. The weak dependence of the 
results on T,^^ means that allowing a time variation of 
T,y^ would not significantly affect our conclusions. The 
value of miimit is sensitive to the timescale over which 
the luminosity decreases. If a very small value for the 
timescale is assumed, as is sometimes the case in the lit- 
erature, then one can obtain apparent sensitivity to a 
very small i^r mass. However, such short timescales are 
unreasonable, given the observed ~ 10 s duration of the 
SN 1987A pulse. 

The supernova parameters are not yet well-known. 
However, the sensitivity with which the neutrino prop- 
erties can be determined using the data from a future 
Galactic supernova depends upon the supernova param- 
eters assumed. We have shown how the supernova and 
neutrino physics can be separately studied using the same 
data, so that the extracted neutrino parameters can be 
determined in an almost model-independent way. 

We assumed that only Vr is massive (note that this can- 
not be distinguished from the case that only is mas- 
sive). If no delay is seen, and nothing further is known, 
then the mass limit obtained would in fact apply to both 
the and Vr masses. If both are massive, then because 
of the quadratic dependence of the delay on the mass, 
the one-mass case is recovered with the larger mass un- 
less the masses are similar. If both are massive and have 
the same mass, the delay would be increased by a fac- 
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tor 2, while the error would be unchanged. This would 
decrease the mass limit by a factor \f2 to about 20 eV. 

If the Vt mass is very large, then the pulse will be so 
delayed and broadened that it will eventually disappear 
below the time-independent background (which for this 
purpose includes solar neutrinos). Assuming that the 
background rate with a lowered threshold is about 10"'^ 
s~^ (about 10 times the solar neutrino rate), the maximal 
mass which can be seen in SNO is about m '--^ 10 keV, 
similar to SK Above this mass, the neutrino would 
be likely to decay over the distance to the supernova, as 
pointed out in Ref. poj . 

Our previous paper [p| showed that using the neutral- 
current excitation of ^"O in SK gives a mass sensitiv- 
ity down to about 45 eV. In SK, one can also use the 
neutrino-electron scattering data, and as noted above, 
that should have a sensitivity down to about 50 eV. If no 
delay is seen, then SNO will place the best limit, about 
30 eV, using neutral-current events on deuterons and cal- 
ibrating the Vf^ data from SK. In that case, the two lim- 
its from SK using neutral-current excitation of ^^O and 
neutrino-electron scattering can only confirm the result. 
For those scattering rates, a mass of 30 eV is insignifi- 
cant, and may as well be zero, so that there is nothing 
to be gained by combining those Signals with that from 
neutrino-deuteron breakup at SNO. On the other hand, 
if a significant delay is seen, then there will be three in- 
dependent determinations of the allowed mass range. 

In conclusion: We have presented a general method, 
including a thorough statistical analysis, of extracting in- 
formation about the possible Vr and i^^ masses from the 
future detection of a Galactic supernova neutrino burst 
by the Sudbury Neutrino Observatory. When such an 
event in fact occurs, the existing mass limits will be vastly 
improved and will approach, or cross over, the cosmolog- 
ical bound. 
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